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Abstract 

We show that the extensively studied equation of state (EOS) of strongly cou- 
pled QED plasma fits the recent lattice EOS data of gluon plasma remarkably 
well, with appropriate modifications to take account of color degrees of free- 
dom and running coupling constant. Hence we conclude that the quark gluon 
plasma near the critical temperature is a strongly coupled color-Coulombic 
plasma. 
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Quark gluon plasma (QGP), the deconfined state of quarks and gluons, is a prediction of 
quantum chromodynamics (QCD). In search of QGP, experiments are on at CERN, BNL etc. 
where heavy ions are accelerated to relativistic energies and made to collide. At sufficiently 
large energy (> few GeV), the fireball formed by collision may be QGP which then expands, 
cools and freezes out into hadrons. These hadrons, photons, leptons etc. are detected and 
analysed to see whether QGP is formed or not. The expansion of QGP, generally described 
by relativistic hydrodynamics, affects the experimental observations. In the description of 
hydrodynamics, the EOS of QGP is needed to complete the set of fluid equations. At present 
all analysis of experimental results are based on ideal equation of state where quarks and 
gluons are noninteracting and pressure or energy density are propotional to the 4 th power 
of temperature. In the theoretical calculations of various signatures of QGP, it is assumed 
that QGP is an ideal plasma. 

However, the recent lattice simulation results show that gluon plasma is not ideal and is 
nonideal even upto five times the critical temperature (T c ) |1J. There were lot of attempts 
to explain this numerical experiment based on the properties of QCD, confinement and 
asymptotic freedom, properties of plasma, etc. All of them, so far, not able to explain the 
data satisfactorily. 

In Ref. ||], they tried to fit the lattice results assuming Coulombic interaction with 
running coupling constant, confinement in the form of bag constant and low momentum 
cut-off in the evaluation of partition function. In Ref. ||, gluons are assumed as quasi- 
particles with mass propotional to plasma frequency. Here also a modified, temperature 
dependent, running coupling constant is used to fit the old lattice results of Ref. ||. As we 
will see it does not fit the recent more refined data [|IJ. Earlier we [|J had assumed Cornel 
potential interaction ||, Coulomb + confinement, between quark and antiquark and used 
Mayer's cluster expansion method to derive EOS. This EOS was used to fit lattice results 
on gluon plasma with partial success. Especially near the critical temperature fitting was 
not good. We suspected that it might be due to the fact that our theory was for weakly 
interacting system and hence valid only in the high temperature limit. To explain the lattice 
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results near T c we need a theory of strongly interacting system. In other words, near T c QGP 
may be a strongly coupled plasma. This is exactly what we find in this letter, namely, QGP 
is a strongly coupled color Coulombic plasma. 

Generally when we say plasma, it is a quasi neutral system of charged particles which 
exhibits collective effects. The so called plasma parameter T, which is the ratio of average 
potential energy to average kinetic energy of particles, is assumed to be weak (<< 1). Lot of 
studies in plasma such as plasma waves, instabilities and other collective effects are in this 
range of T. For V close to 1 and above, it is a strongly coupled plasma (SCP) which modifies 
various properties of plasma J/J. We see that QGP is also a strongly coupled plasma near 
T c . Partially analytic and partially numerical calculations of SCP do exist in the literature 
0. In particular, EOS of SCP is extensively studied and parameterized as a function of T. 

Let us now discuss about QGP. We take QGP as a deconfined, quasi-color-neutral system 
of quarks, antiquarks and gluons with color Coulombic mutual interactions. It is similar to 
QED plasma apart from few modifications due to color degrees of freedom. It is characterized 
by plasma parameter, 

^ < PE > 

<KE>~ T ' { } 

where we have taken the well known Coulombic interactions used in hadron spectroscopy ||. 
The typical value of a s pa 0.5, r av pa lfm and near the critical temperature, T c pa 200 MeV, 
we estimate r ~ 2/3. Hence QGP is a strongly coupled plasma. Compared to QED, the 
fine structure constant, a, is replaced by 4 a s /3 in Coulomb interaction term. r av may be 
estimated as r av = (3/47rn) 1 / 3 and hence 

'47rn\ 1/3 4 a, 
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where 'n' is the number density. Taking n pa ajT 3 we get 



/47raA 1/3 4 , . 

r^Hr 1 o"s, (3) 
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cif is a constant which depends on degrees of freedom. As we discussed earlier there exists 
EOS for strongly coupled Coulombic system as a function of V. Since QGP is also a SCP with 



r given by Eq. (|j), we modify QED nonrelativistic energy density Eqed = (3/2+w ex (T)) nT 
for strongly coupled relativistic QGP (SCQGP) as 



e = (3 + u ex (T))nT , (4) 

where 

_ ^(r) + 3xiQ3r 5 - 7 ^(r) 
} ~ i + 3xio3r" • w 

Here we have taken same u^ e and u^ p as used in SCP and are given by 



uf e (r) = _V3 r 3/ 2 _ 3r 3 



2 ln(3r) + 

8 V 7 2 3 



(6) 
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u ocp = -0.898004 T + 0.96786 T 1/4 + 0.220703 T~ l/i - 0.86097 . (7) 
was derived by Abe || exactly in the giant cluster-expansion theory. 7 = 0.57721... 



is Euler's constant. u^ p was evaluated by computer simulation of one component plasma 
(OCP), where a single species of charged particles embedded in a uniform background of 
neutralizing charges. u ex (T) is derived for strongly coupled Coulombic plasma and we believe 
that it should be valid for any Coulombic plasma except a change, a — > 4a s /3 in T, for color 
Coulombic plasma. Finally we get, 

e(r) = - = 1 + ±u ex (T) , (8) 

where e s = 3a fT 4 , energy density of ideal QGP gas. Using the relation 

dP 

c = T--P, (9) 

we get the pressure 

* T ^¥r{^ +3 £ dTT2e{T) ) /T3 ' (10) 

where P s = dfT 4 and P c is the pressure at critical temperature T c . Temperature depen- 
dence of e(r) seen in lattice data can only come from a s in our model by using running 



coupling constant a s (T). Let us use the running coupling constant (T) used in the lattice 
calculation pi which may be obtained from the modified second order scaling relation, 



Al 



(87r 2 /3/33) 51/121 exp(-47r 2 /V33) , (11) 



N T TX((3) 

where (3 = Q/g 2 and Al,N t are lattice parameters which may be replaced in terms of T c in 
above equation to get, 

£ = (/V/3) 51/121 exp(-47r 2 (/3 c - /3)/33) , (12) 

-L c 

where we have taken A(/3) ~ A(/3 C ), for simplicity. Also it follows from Ref. |T[ that X(f3) 
varies slowly as a function of (3. After some algebra, we get, 

11 (ln(t/t ) + f| ln(127r/33af(t))) 

where t = T/T c and to is 

to = (8vr 2 /3 c /33) 51 / 121 exp(-4vr 2 /? c /33) . (14) 

(t) for given t is obtained by solving Eq. ([il|). t or /3 C is only one parameter and adjusted 
to get the best fitting to lattice data Q. In Fig. 1, we plotted e(t), Eq. (§) and p(t), 

p(t) = y = fa + 3 / ^ 2 e(r)) A 3 , (15) 

and obtained a remarkable good fit to lattice data [I[ for t = .5763 or j3 c = .57. We have 
used a/ = 16 7T 2 /90 which is for gluons. However, the value of (3 C is an order of magnitude 
smaller than that of lattice results. For small value of a s , a^(t) may be approximated as 

a ' {t) ^nEm[ ^ Ht/to) ) ' (16) 

or 

r / x 27T 

«• W - TT« • (17) 

The last expression is similar to the running coupling constant used in Ref. || where to = 
A/T c and also in Ref. |J where to = A/7rT c . A is the QCD scale parameter. For A = 
200 MeV we get T c = 347 MeV and T c = 110 MeV respectively. We have also plotted 



< 18 > 

a measure of nonidealness, and (? s . A reasonable good fit to lattice data is obtained, as 
shown in Fig. 2, except very close to T c . In Fig. 3, we plotted (t) and T(t) as a function 
of t. We see that even at T = 5T C , gluon plasma is strongly coupled plasma and running 
coupling constant is not small. 

It is interesting to compare earlier few theories with our present theory. In Fig. 4, the 
best fits for p(t) for three models (Ref. |5J, Ref. [f|] and present model) are plotted. As we 
see the fitting of lattice data with the theory improves as we go from quasi-particle theory, 
our earlier result and present result. 

In conclusion, we obtained a remarkably good fit to the recent lattice result [p]] using 
EOS of strongly coupled plasma of QED with a modification for color degrees of freedom 
and with a running coupling constant which was used in lattice calculations. Hence we 
conclude that QGP is a strongly coupled Coulombic plasma of color charges. Just like in 
SCP, SCQGP may have a dramatic effects on various signatures, collective effects and other 
properties. They need to be recalculated before we make confirmation of QGP in relativistic 
heavy ion collisions. For e.g., screening length decreases rapidly near critical temperature 
and we may have serious consequences in J/ip suppression results. 

I thank P. K. Kaw and J. C. Parikh for useful discussions on strongly coupled plasma. 
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FIGURES 

FIG. 1. Plots of p(t) = P/P s (dashed curve) and e(t) = e/e s (continuous curve) as a function 
of t = T/T c from our model and lattice results (dots), respectively. 

FIG. 2. Plots of A = (e — 3P)/e s (continuous curve) and c 2 s (dashed curve) as a function of 
t = T/T C from our model and lattice results (dots), respectively. 

FIG. 3. Plots of ag (t) (dashed curve) and T(t) (continuous curve) as a function of t. 

FIG. 4. Plots of the best fits for p(t) from three models, Ref. [3] (dotted curve), Ref. [4] 
(dashed curve), present model (continuous curve) and lattice results (dots) as a function of t. 
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